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Search for the Electric Dipole Moment of the r Lepton 

K.Inami for the Belle collaboration 

'^Department of Physics, Nagoya University, 
Furo-cho, Chikusa-ku, Nagoya, 464-8602, Japan 

We have searched for a T/CP violation signature arising from an electric dipole form factor {dr) of the r lepton 
in the e'^e~ t^t~ reaction. Using an optimal observable method for 29.5 fb~^ of data collected with the Belle 
detector at the KEKB collider at ^/s — 10.58 GeV, we obtained the preliminary result Re^dr) = (1.15 ± 1.70) x 
10~^^ecm and Im{dr) = (—0.83 ± 0.86) x 10~^^ecm for the real and imaginary parts of dr, respectively, and set 
the 95% confidence level limits -2.2 < Re{dr) < 4.5 (10""ecm) and -2.5 < Jm(dr) < 0.8 (10""ecm). 
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While large CP violating asymmetry in the 
quark sector has recently been confirmed in B- 
meson decay at B factories the Standard 

Model (SM) does not predict any appreciable 
CP violation (CPV) in the lepton sector. How- 
ever, physics beyond the SM, such as multi-Higgs, 
SUSY, and lepto-quark models |^,|| , could result 
in CPV in leptonic processes; since new bosons 
and Higgs would strongly couple with heavy par- 
ticles through quantum loop effects, CP violating 
phenomena could be significantly enhanced, par- 
ticularly in r decay, due to its very large mass 
compared to other leptons. 

The contribution from CP violating interac- 
tions in the r-pair production process can be 
parametrized in a model-independent way using 
the electric dipole moment (or form factor) for 
the T lepton. Up to now, searches for the electric 
dipole form factor dr have been performed at 
LEP 1^,^ in the reaction e"''e^ t^t^^ and 
by ARGUS in the e+e" t+t" reaction 
and the following upper limits at 95% confidence 
level were obtained: —3.1 < dr < 3.1 (lO^^^ecm) 
from L3 g, -3.8 < rf^ < 3.6 (IQ-i^ecm) 
from OPAL § and \Re{dr)\ < 4.6 x lO-^^ecm, 
\Imidr)\ < 1.8 x lO-i^ecm from ARGUS Q. 

We have searched for a form factor dr at the 
jTT vertex in the process e+e^ j* ^ t^t^ us- 
ing triple momentum and spin correlation observ- 
ables. The effective Lagrangian can be expressed 



Ccp = -idri-sya'^" j5Tdf,A^ 



(1) 



where the electric dipole form factor dr depends 
in general on s, the squared energy of the t- 
pair system. In common with other analyses 
we ignore this possible s-dependence, assuming 
dr{s) = dr, constant. The squared spin density 
matrix (A^p^od) ^o'' the reaction e+(p)e^(— p) 



T+ik,S+)T-{-k,S-) is given by [§ 



M 



prod 



= MIm + Re{dr)M%, + Im{dr)Mj„ 



(2) 



MIm = 

4 

[fco + ml + Ifc^l(fcp)^ - S+S-\kf{l - (kpf) 

fcf, 

+ 2{kS+){kS-){\kf + (fco - mrf(kpf) 



- 2fco(fco - mr)(kp){(kS+){pS-) + {kS-){pS+)) 

+ 2kl{pS+){pS^), 

3 

(rriT -I- (fco 



Ml, = ^^\k 
fco 



mr){kpf){S^ 



+ ko{kp){S+ 
fco 

+ ko{kp){S^ 



X S-)p], 
[—{rrir + (fco 
-S-)p], 



mr){kpf)iS, 



(3) 
X S-)ie 

(4) 
- 5-)fe 

(5) 



where fco is the energy of the r, is the t mass, 
p is the momentum vector of e^, k is the momen- 
tum vector of r"*" in the center-of-mass frame, S± 
are the spin vectors for , and the hat denotes a 
unit momentum. We disregard the higher order 
terms proportional to jd^j assuming dr{s) to be 
small. TWg,^ corresponds to the SM term. M%^^ 
and A^/to are the interference terms (related to 
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the real and imaginary parts of dr, respectively) 
between the SM and CPV amplitudes. is 
CP odd and T odd, while Mj^ is CP odd, but 
T even. In the above equations, e+ and e~ are 
assumed to be unpolarized and massless particles. 

We adapt the so-called optimal observable 
method |^ , which maximizes the sensitivity to dr 
by optimizing the relevant quantities, O^e and 
Oim, for measured kinematic variables, as de- 
scribed below. Here the optimal observables are 
defined as 



O 



Re 



•ML 



(6) 



The mean value of the observables is expressed as 



(ORe) CX / OReda CX / OReM 



'^prod 



: J MlJ(l) + Re{dr) J 



(7) 

where the integration is over the phase space 
((/)) sustained by the relevant kinematic variables. 
The cross-term containing the integral of the 
product of ^^^'^ -M/m drops out because of 

their different symmetry properties. The expres- 
sion for the imaginary part is similar. The means 
of the observables {Ore) and {Ojm) are thus ex- 
pressed as linear functions of dr 

(Orc) = a_Re • Re{dr) + bRe, 

(Olm) = a/m • Im{dr) + bl„i. (8) 

Eight different final states in the decay 
of T-pairs, {evv){^i>v), {evv){-nv), {^i>v){-kv), 
{evv){pv), {iivD){pv), {■Kv)[pv), {pv){pv), and 
{■Kv){'Kv), are analyzed, where all particles except 
V and v are positively or negatively charged. Be- 
cause of the undetectable particles, we can not 
fully reconstruct the quantities k and S± . There- 
fore, for each event we calculate possible kine- 
matic configurations and obtain the mean value 
of -^sM' -M/je ^"^d WI/to by averaging over the 
calculated configurations. In the case when both 
r leptons decay hadronically (t ttv or pv) , the 
r flight direction is calculated with a two-fold am- 
biguity and we take the average of Algj^j, Mr^ 
and A^/^ over the two solutions. In the case 
when either one or both r leptons decay lepton- 
ically (r — > ei'D or pvv), a Monte Carlo (MC) 



treatment is adopted to take into account the ad- 
ditional ambiguity in the effective mass of the vv 
system {m,^p). For each event we generate 100 
MC configurations using a hit-and-miss approach 
by varying m^p, and compute the averaged Alg^, 
A^^g and over successful tries in which the 

T direction can be constructed kinematically. In 
the calculation, we ignore the effect of undetected 
photons coming from initial state radiation, ra- 
diative r decays and bremsstrahlung. This ef- 
fect has been examined and found to be negligible 
compared to the total errors. 



2. Data and event selection 

In this analysis, we used 26.8 million r-pairs 
(29.5 fb~^) accumulated with the Belle detec- 
tor [|o| at the KEKB accelerator KEKB 
is an asymmetric energy e'^e" collider with a 
beam crossing angle of 22 mrad. Its center- 
of-mass energy is 10.58 GeV, corresponding to 
the T(4S') resonance, with beam energies of 8 
and 3.5 GeV for electrons and positrons, respec- 
tively. Belle is a general-purpose detector with 
an asymmetric structure along the beam direc- 
tion. Among the detector elements, the central 
drift chamber (CDC) and the silicon vertex detec- 
tor (SVD) are essential to obtain the momentum 
vectors of charged particles. The combined infor- 
mation from the silica Aerogel Cherenkov coun- 
ters (ACC), the time-of- flight counters (TOE), 
the Csl electromagnetic calorimeter (ECL), and 
the /i/KL detector (KLM) is used for particle 
identification. 

The MC event generators KORALB / 



TAUOLA 1 12 are used for r-pair production and 
decays. The detector simulation is performed by 
a GEANT-based program, GSIM. Actual data 
and MC generated events are reconstructed by 
the same program written by the Belle collabo- 
ration. 

We chose the eight final state r-pair modes 
mentioned above. All of the final state particles 
were reconstructed with the following conditions. 
Each charged track is required to have a trans- 
verse momentum pt > 0.1 GcV/c. Photon candi- 
dates should deposit an energy of £' > 0.1 GeV 
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Table 1 ^ 

Yield, purity and background rate obtained by the event selection described in the text, where the purity 
was evaluated by MC simulation and its error comes from the MC statistics. 





Yield 


Purity (%) 


Background mode (%) 




250,948 


96.6 ±0.1 


27- 


pp{1.9), TT e7r(l.l). 


en 


132,574 


82.5 ±0.1 


TT - 


^ ep(6.0), eX(5.4), e^(3.1), eX*(1.3). 




123,520 


80.6 ±0.1 


TT - 


^ pp{5.7), pK{5.3), pp{2.%), 27 pp{2.Q). 


ep 


240,501 


92.4 ±0.1 


TT ~ 


-> e7r7r07r"(4.4), eK*{l.7). 


MP 


217,156 


91.6 ± 0.1 


TT - 


-> ^7r7r07rO(4.2), pK*{l.&), 7rp(1.0). 


Tip 


110,414 


77.7 ±0.1 


TT - 


^ PP(5.1), Kpi'i.Q), TTTTTT^TT^i.S), pp{2.7). 


PP 


93,016 


86.2 ±0.1 


TT - 


^ P7r7r07r0(8.0), pK*{'i.\). 


vrvr 


28,348 


70.0 ±0.2 


TT - 


-y 7rp(9.2), ttK{9.2), 7r^(4.7), 7ri^*(2.0). 



in the ECL. A signal event is required to have 
two charged tracks with zero net-charge and no 
photon apart from p^ tt^tt", tt^ — > 77. 

A track was identified as an electron using a 
likelihood ratio combining dE/dx in the CDC, the 
ratio of energy deposited in the ECL and momen- 
tum measured in the CDC, the shower shape of 
the ECL and the hit pattern from the ACC. The 
identification efficiency was estimated to be 92% 
with a TT* fake rate of 0.3% jl^. A muon was 
evaluated from its range and hit pattern at the 
KLM detector; its efficiency and fake rate were 
obtained by MC as 91% and 2%, respectively. A 
pion track was found by requiring that a track 
be identified as a hadron by the KLM informa- 
tion, and not identified as an electron. The ef- 
ficiency was estimated to be 81%. The purity 
for the obtained samples is about 89%. A p^ 
was reconstructed from a charged track and a tt'^ 
where the track should be neither an electron nor 
a muon, and for tt" 77, the reconstructed 7r° 
should have an invariant mass between 110 and 
150 MeV/c^ and a momentum in the laboratory 
frame larger than 0.2 GeV/c. In order to sup- 
press the background and to enhance the par- 
ticle identification ability, the detection of lep- 
tonic particles was restricted to within the bar- 
rel region, -0.60 < cos 6'^'^'^ < 0.83, and that 
of single pions to within the KLM barrel region, 
-0.50 < cose^^^ < 0.62, where O^''^ is the po- 
lar angle relative to the e~ beam direction. For 
the same reason, the particle momentum was re- 
quired to be greater than 0.5 GeV/c for an elec- 
tron, 1.2 GeV/c for both a muon and pion, and 



1.0 GeV/c for p^ in the laboratory frame. 

The dominant backgrounds are due to two- 
photon as well as Bhabha and pp processes. In 
order to remove two-photon events, we required 
the missing momentum not to be directed to- 
wards the beam-pipe region (imposing a selec- 
tion -0.950 < cosO^'^^ < 0.985), and to reject 
the latter processes we required that the sum of 
the charged track momenta be less than 9 GeV/c 
in the center-of-mass frame. Additional selections 
were imposed particularly on the ctt mode where 
a large number of Bhabha events could contribute 
through misidentification. For the ctt mode, we 
remove events which satisfy the following crite- 
ria: the opening angle of the two tracks in the 
plane perpendicular to the beam axis is greater 
than 175°, and their momentum sum is greater 
than 6 GeV/c in the r-pair rest frame. Finally, 
we removed events in which the r flight direction 
could not be kinematically reconstructed, which 
mostly arise from r-pairs having hard initial-state 
radiation and misidentified r-pair backgrounds. 

The yield of events passing this selection is 
given in Table |] for each of the eight selected 
modes. The mean energy of the r-pair system in 
the obtained sample is ^/s = 10.38 GeV; this sets 
the scale at which dr{s) is measured. Because 
of events with soft radiated photons, the energy 
scale is slightly lower than the beam energy. The 
dominant background sources are also listed in 
the table. Hadronic r decays with two or more 
tt'^'s make a large contribution of a few percent. 
For the modes including tt^ , the misidentification 
of kaons and muons as pions yields other back- 
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Figure 1. Momentum distributions of (a) e^, (b) 
fj,^ , (c) TT^, and (d) in the laboratory frame. 
The points with error bars are the data and the 
histogram is the MC expectation. The latter is 
scaled to the total number of entries. The hatched 
histogram is the background distribution evalu- 
ated by MC. 



grounds; for example, 5.3% of fiK final states 
is included in the fiir mode. The other back- 
grounds are estimated by MC to be a few percent 
from two-photon processes, and less than 1% from 
Bhabha, /i/i, and multihadronic processes. 

Figs and show the resulting momen- 
tum and cos 61' distributions, respectively, for 
charged particles in the laboratory frame. Very 
good agreement with MC is found, except for 
low-momentum electrons (Fig. |l|(a)) and pious 
(Fig. 0(c)). The dip in the cos6' distribution of 
the muon (Fig. ||(b)) is due to an efficiency drop 
at the overlapping region between the barrel and 
endcap KLM elements. 

The resultant Ors and O/m distributions are 
shown in Fig. ^ along with those obtained from 
MC simulation with dr = 0. Good agreement is 
found between the experimental data and the MC 
samples. The ratio of the data to MC shows flat 




Figure 2. The cos 6* distributions of (a) e^, (b) 
/i^, (c) 77^, and (d) in the laboratory frame. 
The meanings of the points and histograms are 
the same as in Fig. |[ 



distributions around 1.0, although they are not 
shown here. 

3. Extraction of dr 

In order to extract the d-r value from the ob- 
servable using Eq. (^), we have to know the co- 
efficient a and the offset b. In the ARGUS anal- 
ysis [0 which also used the optimal observable 
method, the first term of Eq. (|^) was assumed to 
be negligible because of the property of the CP 
violating term under the ideal detection hypoth- 
esis J Aij^^dcj) = 0. The value of dr was thus 
obtained as the ratio of the observable's mean to 
the second term, Re{dr) = {O Re) / (Oj^^) . Exper- 
imentally, the detector acceptance r/ affects the 
means of the observables as 

{One) cx J r,{4>)OReM%,^dc^. (9) 

Similar expressions are obtained for the imagi- 
nary part. This means that the first term of 
Eq. (1^) is not necessarily zero and the coefficient 
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Figure 4. Re{dr){I'm{dT)) dependence of the mean 
of the observable {ORe){{Oim}) for the irp mode. The 
closed circles show the dependence for Re{d-r) and the 
open circles show the dependence for Imld-r). The 
lines show the fitted linear functions. 



Figure 5. Sensitivity a (a) and offset h (b) for each 
mode from MC. The closed circles show the parame- 
ters for Re{dT) and the open circles show the param- 
eters for Im(dT). The errors are due to MC statistics. 



may differ from (C|jg) when the detector accep- 
tance is taken into account. In their study, the 
acceptance effect produced the largest systematic 
uncertainty, on the order of 10~^^ecm. 

In order to reduce this systematic effect, 
we use both parameters a and b extracted 
from the correlation between {ORe){{Oim)) and 
Re{dr){Im{dr)) obtained by a full MC including 
acceptance effects. An example of the correla- 
tion between (C'fle)((C/m)) and Re{dr){Im{dr)) 
is shown in Fig. ^. Each point is obtained from 
MC with detector simulation and event selection. 
By fitting the correlation plot with Eq. (|8|), the 
parameters a and h are obtained. The misiden- 
tified background from r-pair events shows some 
dependence on dr, because the spin direction is 
correlated with the momenta of the final state 
particles. Therefore, both the coefficient a and 
offset h are corrected with the parameters ob- 
tained from the MC misidentified background. 
The resultant coefficients and offsets are shown 
in Fig. |. 



As can be seen from the values of the coeffi- 
cient a showing the sensitivity to d,-, the irp and 
PP modes have the highest sensitivities, while the 
TTTT mode yields a somewhat lower value. This is 
because both S and k in the tttt mode are inferred 
only from the tt momentum vector, so that tak- 
ing the mean value of two possible solutions with 
equal weights significantly reduces the achievable 
sensitivity. On the other hand, the momentum 
sum and angular distribution oi p tttt^ decay 
in the -np and pp modes provide different weights 
on the possible solutions, and therefore there is 
less reduction in sensitivity. Other modes, which 
include leptons, have low sensitivities because of 
a lack of information about the directions and 
spin orientations due to the additional missing 
neutrinos. 

Non-zero offsets him are seen for the imaginary 
part, due to the forward/backward asymmetry in 
the acceptance of the detector. 
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Table 2 

Systematic errors for Re{dr) and Im{dr) in units of 10~"'^^ecm. 



Re[dr) 






flTT 


GO 


HO 


7TO 


oo 


TTTT 


Mismatch of distribution 


0.80 


0.58 


0.70 


0.11 


0.15 


0.21 


0.16 


0.06 


Charge asymmetry 


0.00 


0.01 


0.01 


0.01 


0.01 


0.01 






Background variation 


0.43 


0.12 


0.07 


0.07 


0.08 


0.03 


0.04 


0.05 


Momentum reconstruction 


0.16 


0.09 


0.24 


0.04 


0.06 


0.06 


0.04 


0.45 


Detector ahgnment 


0.02 


0.02 


0.01 


0.00 


0.01 


0.01 


0.02 


0.03 


Total 


0.92 


0.60 


0.74 


0.14 


0.18 


0.22 


0.17 


0.46 


Im{dT) 


e/i 




flTT 


ep 


pp 


TTp 


PP 


TTTT 


Mismatch of distribution 


0.43 


0.02 


0.05 


0.12 


0.01 


0.05 


0.10 


0.41 


Charge asymmetry 


0.13 


0.44 


0.43 


0.02 


0.09 


0.15 






Background variation 


0.08 


0.07 


0.02 


0.01 


0.03 


0.02 


0.03 


0.06 


Momentum reconstruction 


0.03 


0.03 


0.06 


0.00 


0.02 


0.02 


0.04 


0.04 


Detector alignment 


0.01 


0.03 


0.02 


0.01 


0.01 


0.02 


0.01 


0.05 


Total 


0.46 


0.45 


0.44 


0.13 


0.10 


0.16 


0.11 


0.42 



4. Systematic uncertainties 

Although in general the MC simulation repro- 
duces the observed kinematic distributions well, 
the disagreement seen in Figs. |l| and ^ causes 
an error dominating the systematic uncertainty. 
Its uncertainty with respect to d^ was studied 
by reweighting the MC distributions by the ra- 
tio of data to MC. The second significant un- 
certainty originates from possible charge asym- 
metry in the detection efficiency. The ratio of 
yields, N{a+ P') /N[a- (}+), for data and MC 
agrees within 1% accuracy, where a and (3 are the 
relevant charged particles from the t decays. This 
systematic uncertainty was evaluated by varying 
the detection efficiency by ±1%. The effect is of 
the same size as the statistical error for Im(dr), 
while it is negligible for Re{dr). The backgrounds 
provide additional systematic errors for dr be- 
cause the parameters a and h are corrected for 
the background contributions, which are assessed 
by varying the background rate. The influence 
of the momentum reconstruction of charged par- 
ticles and photons was also checked by applying 
scaling factors corresponding to the difference be- 
tween the data and MC. In order to examine a 
possible asymmetry arising from the alignment 
of the tracking devices, we measured the differ- 
ences of the polar angles, AO = 9^^^ — 0^^, and 
the azimuthal angles, A(p = — ^^M' between 



two tracks in e+e^ P^ P-^ events, and found 
a small deviation from a back-to-back topological 
alignment in each direction as AO — 1.48 mrad 
and A<f> = 0.36 mrad. By applying an artificial 
angular deviation of this magnitude to one of the 
charged tracks, the residual value of the observ- 
ables is calculated to be negligible compared with 
the other errors. The systematic errors are listed 
in Table |. 

5. Result 

The values of dr extracted using Eq. (^) are 
listed in Table ^ along with the statistical and 
systematic errors. All results are consistent with 
zero within the errors. 

Finally, we obtain mean values for Re{dr) 
and Im{dr) over the eight different t+t" modes 
weighted by quadratically summed statistical and 
systematic errors. The resultant preliminary re- 
sult for the electric dipole form factors are 

Re{dr) = (1.15 ± 1.70) X lO^^ecm, (10) 
Im{dr) = (-0.83 ± 0.86) X lO^i^ecm, (11) 

with corresponding 95% confidence limits 

- 2.2 < Re{dr) < 4.5 (lO^^ecm), (12) 
-2.5 < Im{dr) < 0.8 (10-i^ecm). (13) 

This investigation has improved the sensitivity 
to the r lepton's electric dipole form factor by 
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Table 3 

Results for the electric dipole form factor. The 



nrst 
atic. 


error is statistical and 


the second is system- 


1. 




Re{dr) (10-i'^ecm) 


Im{dr) (10-i^ecm) 


2. 


e/i 


2.25 ± 1.26 ±0.92 


-0.41 ±0.22 ±0.46 


eir 


0.43 ±0.64 ±0.60 


-0.22 ±0.19 ±0.45 


3. 


fXTT 


-0.41 ±0.87 ±0.74 


0.15 ±0.19 ±0.44 


ep 


0.00 ±0.36 ±0.14 


-0.01 ±0.14 ±0.13 


4. 


HP 


0.04 ±0.42 ±0.18 


-0.02 ±0.14± 0.10 


TTp 


0.34 ±0.25 ±0.22 


-0.22 ±0.13 ±0.16 




PP 


-0.08 ±0.25 ±0.17 


-0.12 ±0.14 ±0.11 




TTTT 


0.42 ± 1.17 ±0.46 


0.24 ±0.34 ±0.42 


5. 


All 


0.115± 0.170 


-0.083 ±0.086 





a factor of more than 10 compared to previous 
measurements. 
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Figure 3. Distributions of the optimal observables Ore and Oim for each mode. The upper figure of each mode 
is for Oije and lower figure is for Oirn- The closed circles are the experimental data and the histogram is the 
MC expectation with cij- = 0, normalized to the number of entries. The hatched histogram is the background 
distribution evaluated by MC. 



